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Abstract 

We imitate the spectrum character of one-dimensional disorder system with our new rain model. 
It has been shown that the transmission spectrum can be approximately characterized by the 
model, which include some coupled lorentzian transmission peaks on the spectrum and embody 
the coupling of the local modes in the disorder system. 
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The transport of optical waves through disordered systems can exhibit remarkable inter- 
ference effects, in analogy with the transport of electrons in solids. This leads to interesting 
optical phenomena[l| of which the most surprising is that of Anderson localization of light j^J. 
The observation of this extraordinary effect in three dimensional disordered optical systems 
requires very strong scattering that can be achieved only in selected materials [3j]. For lower 
dimensional systems, however, the situation is different. In one and two dimensional dis- 
ordered systems, localization can always be reached for a sufficiently large sample size J4] 

ly localized in ID random systems, even though 
6[ and Tartakovskii [7] predicted that even in one 
dimensional (ID) localized systems nonlocalized modes exist that extend over the sample via 
multiple resonances. These nonlocalized modes, called necklace states(NS), have a transmis- 
sion coefficient close to 1 and become extremely rare upon increasing the sample thickness. 
Nevertheless, they dominate the average transmission coefficient, even at large thickness. 
Evidence for the existence of necklace states was recently found in time-resolved trans- 
mission experiments js], and subsequently also observed in experiments with microwaves^. 
It is noted that these works about NS are mainly concentrated on the one dimensional 
systems for its rather easily calculation and manufacture. The basic method and ways is 



However, not all modes are exponential 
Anderson localization occurs. Pendry 5). 



commonly to analyze the transmission spectra of the system whether experimentally lOl. |ll| 



or theoretically 
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since the transmission spectra is essential for ID disorder system 
and reveal important formation of the system. For example,when the phase of the transmis- 
sion coefficient is the integral multiple of 2ir, it corresponds to the order of NS 121] and the 
platforms in the spectrum usually can be used to identify the occurrence of NS. Since the 
transmission spectra directly implicate the intrinsic properties of the system, It is reasonable 
that one can artificially construct transmission spectra according to some aforehand rules 
inferred form the physical evaluation. Some statical properties of physical variables of the 
system can be explored based on the constructed spectra. It is equivalent that one is in- 
quiring interactively the relations between the statical properties and the possible physical 
mechanism. Based on the above idea, In this letter we present a model called after rain 
model to imitate the transmission spectrum of the one dimensional disorder system. Gen- 
erally, the transmission spectra of localized ID systems exhibit many randomly distributed 
high-transmission peaks. These high-transmission peaks originate from resonances of the 
system with localized modes Ijj] and result in big fluctuations in the transmission coefficient 
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FIG. 1: schematic figure of our model 
t[i[ [isj] . The localized modes decay exponentially and the ensemble average of over many 



realizations of the disorder decays linearly with the sample thickness L jo]. 

Our model is as follows. Consider a one-dimensional disorder system with length L. The 
localized states inhabiting in the system corresponding to the eigenfrequency u n usually have 
the form E n oc e~\ x ~ Xn \l^ except a oscillating factor, where £ is the localization length and 
x n . For simplicity, omitting the fine structure of the field for the moment we assume the field 
localized at x n has simply the form of the envelope. Since the necklace states are usually 
constructed by a series of localized states close to each other spatially so that they seem 
like the necklace and the resultant transmission spectrum manifest a high lorentzian peaks 
bundle. Therefor, first of all we should investigate what the transmission spectrum would be 
like if the fields in system are superposed by multiple localized states. Here we have to use 
experiential conclusion that the transmission coefficient t is expressed as t = ^a" {g" (o)' g" (L) } 
With the above approximation and assuming x n < L/2, the transmission coefficient for 
a local state x n is t = fr^y = e~ L ^e 2Xn ^. For the field of iV coupled states E N (x) = 

N N 

Yln=i e~\ x ~ Xn \/£, the corresponding transmission coefficient = = S tr/ 2 / Y^, tn 1 ^ 2 - 

n n 

Since the transmittance peaks in the bundle are lorentzian and has the form T(oS) oc 
, _ \2\_1\ / 9 \2 ; where the center frequency u n and all the random variables and 

the halfwidth Au can be determined by quality factor Q with relation Q = Accord- 



ing to the definition of Q = ^ ^ored energy _ w store. _ wq_ ^gi halfwidth can be 

to ^ u p OW6r \ QSS u p loss Aw — lj) 



written as Aw = w l ° aa . In one dimensional system how is the quality factor correlated 
with the transmittance? Since quality factor Q is the measure of the sharpness of re- 
sponse of the cavity to external excitation, it's reasonable that the system is regarded a 
cavity enclosed by the vacuum and the transmitted energy corresponds the power loss. 
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Here consider a imaginary pipe with unit area though the system. So the power loss 
can be expressed transmitted average energy current (S) x= l = while the stored en- 
ergy in the pipe W st0 re = J(U)dv = J E Re(e\E\ 2 )dx. It is found from above results 
that Atu = — E ^\j = cE } L li , where c is the light speed, Z = \/ Unit® is the vac- 

2Z efE 2 dx 2t r J E 2 dx^ to r i v V r-vi u 

uum impedance, e r is the average relative permittivity of the system. For the above field 
E n (x) = e~l x-x "l^, we obtain 

e 4 *"/* 2vrc 

U ~ ~ (e 2L /5 + e 4*„/£ _ 2e 2 ( L+:E ")/«)27re r f U 

To imitate the spectrum of the necklace states we assume many lorentzian peaks randomly 
fall on the frequency axis with random center frequency u n , which is related with another 
spatial random variable x m that is the center of local states. For each (uu n ,x m ), we first 
construct the transmittance coefficient 

e ( 2 ^- L )/ 2 «(Aw/2) 
t nm {u>, u n , x m ) - (2) 
y/{u - Cd n ) 2 + (Aw/2) 2 

and the total transmittance spectrum can be obtained 

V t 1/2 

Hotal{U) - — ~Zy2 \ 6 ) 

To explore the possible mechanism we introduce the repulsive interaction between these 
'frequency raindrops' which are characterized by the corresponding overlap integral between 
the two localized states. Specifically speaking, for the states localized at Xi,Xj, the over- 
lap integral is = J^e^ x ~ Xi ^ x ~ X: >^^dx = (£ + Ax)e~ Ax ^ and the moduls square 
J™ \E(x)\ 2 dx = £. In the rain process the frequency points that satisfy the condition 
Iij < at; are permitted to be coupled into the t tota i, where parameter a control the coupling 
strength. We call the process as two order process since it only involves two coupled fre- 
quencies and two localized states. To check the validity of the operation, we compared the 
correlation curves to frequency interval with those obtained from the real one dimensional 
disorder system by transfer matrix method. The correlation is defined as: 

= (OnT(q; + Au) - (lnT(w + Au)))(\nT(u) - (lnT(u)))) 
y/(QnT(u + Aw) 2 - (lnT(w + Aw)) 2 )(lnT(w) 2 - (lnT(w)) 2 )) 
, where (•) means the average over all realizations of disorder. For the two order process 
we introduce a coupled mode theory to phenomenologically characterize the effect on the 
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FIG. 2: (a)lnT(w) spectrum for Transfer method(solid) and rain model (dashed) (b) Correlation for different frequency 

interval 

transmission coefficient by the fact that mode coupling transform the two random frequencies 
u m ,u n into the renormalized frequencies u' m ,u' n : 

d 2 a± dai 



+ 7i-7T + Wia.i = <?i2«2 + e 



-iw't 



dt 2 dt 
d 2 a 2 da 



+ l2—rr + u 2 a 2 = q 2 \a x 



(5) 



dt 2 1 dt 

Where ai(t) = a w e~ lU!lt , a 2 {t) = a 20 e~ ia;2 * and a 10 , a 2 o, w' is to be determined. According 
to the condition of nonzero solution of a w , a 2 o, 



(6) 



we can get two renormalized peaks frequency u[, u>' 2 with positive real parts, where 
the factor e~ tu ' 1 describes the excitation of incident wave and for simplicity the coupling 
coefficients q\ 2 = 921 = q oc Iij, the damping factor 7„ oc Acu 

In Fig. 2 we compare the lnT(a;) and the correlation obtained by transfer matrix method 
for a 3000 layers disorder system with d\ = 200nm, €a — 1, di — lOOnm, e# = 2 with the 
corresponding results obtained by the rain model. In Fig. 2(b) for one order approximation 
where lorenztian peaks are summed directly the correlation by the rain model is only in 
good agreement with the system with low random strength where w = 0.04. However for 
higher random strength the one order approximation apparently fail in approximating the 
real disorder system. For the real system our calculation show that when random strength 
w grow gradually the correlation for a fixed frequency interval rise rapidly and descends 
slowly after passing by a maximum, which is shown in Fig. 3. Here for our the calculated 
3000 layers' system the critical random strength w c ~ 0.12. It is found that when more and 
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FIG. 3: Correlation curve corresponding to different random strength w for some fixed frequency intervals 

more lorenztian peaks are coupled into the transmission coefficient spectrum, which means 
the higher order processes, the whole correlation curve also ascend rapidly as shown by the 
green dashed dotted line in Fig.2(b). For a given random strength w the correlation of the 
real system can be always approximated by the rain model with specific number of order 
process because it can be raised consecutively with the increment of number of order. For 
example, for w = 0.054 rain model can give a good approximation with ~ 20000 peaks and 
a = 0.5, L/C, = 8 which is marked with blue dashed line in Fig.2(b). The corresponding 
spectrum is shown in Fig. 3(a). It is noted form Fig. 3(a) that in the range 15.97 ~ 15.99 the 
spectrums keep the close similarity and some double peaks structure begin to occur, which 
imply the thrown coupling double peaks play a important role in the similarity of correlation 
curve. 

In conclusion, we present a model to simulate the transmittance spectrum of the one 
dimensional disorder system. The model describe the construction of the spectrum based 
on a thought that many lorenztian peaks can be superposed in a way embodying some 
interaction between these single peaks. Furthermore, we introduce overlap integral between 
localized states and the coupling of modes to characterize the interaction phenomenologically. 
In the above way lorenztian peaks randomly drops on the frequency axis like raindrops and 
for each peak the characteristic parameter is determined by a localized states with a random 
position. Here it involves two aspects of randomness: frequency and spatial domain. During 
the sampling of the random variables uniform distribution is used. In the weak disorder 
the model can give rather better approximation. Because the model associate the form of 
transmittance spectrum with the interaction between single peaks and the position of the 
localized states it provide a intuitionistic tool for the subject that depends on the analysis 
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of the transmittance spectrum such as the study of the necklace states. 
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